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Abstract 

We demonstrate a minimally invasive nuclear magnetic resonance (NMR) 
technique that enables determination of the surface-area-to-volume ratio 
(S/V) of soft porous materials from measurements of the diffusive exchange 
of laser-polarized ^^^Xe between gas in the pore space and ^-^^Xe dissolved in 
the solid phase. We apply this NMR technique to porous polymer samples 
and find approximate agreement with destructive stereological measurements 
of S/V obtained with optical confocal microscopy. Potential applications of 
laser-polarized xenon interphase exchange NMR include measurements of in 
vivo lung function in humans and characterization of gas chromatography 
columns. 
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Porous media are ubiquitous in nature, e.g., granular materials, foams, ceramics, oil- or 
water-bearing "reservoir" rocks, and animal and human lungs. Determining the structure of 
these materials is relevant to a wide range of scientific and technological problems, ranging 
from the coarsening of foams to the transport properties of subsurface fluids to cardiopul- 
monary physiology and medicine. In this paper we demonstrate a noninvasive technique to 
characterize the the surface-area-to-volume ratio {S/V) IQ of "soft" porous media: i.e., ma- 
terials in which there is significant gas solubility in the solid phase. Such a description applies 
to many different materials, including porous polymer granulates used for filtering, porous 
polymer bead packs used for radio-immunoassay, resin columns used for chromatographic 
separation, and lung and sinus tissue in humans and animals. 

S/V influences numerous interactions within porous media, including fluid diffusion and 
transport, electric charge distribution, and chemical exchange. For example, S/V is impor- 
tant in determining pulmonary function since the lung is the site of O2 and CO2 exchange 
between the body and the external environment. There are numerous techniques to mea- 
sure S/V in porous media, including stereology with traditional and confocal microscopy, 
light scattering , and mercury intrusion porosimetry . However, only nuclear magnetic 
resonance (NMR) provides a noninvasive and nondestructive technique for determining 3D 
structure, including S/V 0. 

In particular, NMR measurements of the time-dependent diffusion coefficient of a liquid 
filling the pore space has been used to determine S/V in a variety of porous media [|^,|^ 
including glass beads 0, reservoir rocks [jlO[, and biological samples [|Tl|. Recently, we 



extended this technique to NMR measurements of xenon gas diffusion within porous media 
T^ , p!3[| , thereby probing longer length scales than with traditional liquid-infused diffusion 



measurements [^|10| because of fast gas diffusion and long xenon spin polarization lifetimes. 
Gas diffusion NMR may also be used to study porous media that are not compatible with 
water-saturation, such as non-wetting polymers and in vivo lungs. However, for materials 
with pore size less than about 300 /im (such as lung alveoli), the high diffusivity of gases 
causes systematic errors in the determination of S/V because the rms diffusion distance 



during the application of RF and magnetic field gradient pulses is comparable to the mean 
pore size |T^-|14|. Reducing this effect with faster gradient pulses, buffer gases, or higher 
gas pressures is impractical [|1^ . 



Thus we have devised an alternative NMR method ( "xenon interphase exchange NMR" ) 
to measure S/V in soft (i.e., xenon-soluble) porous media. In this technique we determine 
S/V from the rate of ^^^Xe polarization (not chemical) exchange between gas in the pore 
space and ^^^Xe dissolved in the solid. First, we infuse the sample with laser-polarized (i.e., 
highly spin-polarized or hyperpolarized) ^^^Xe gas. Next, we rapidly destroy ("quench") 
all dissolved-phase ^^^Xe polarization using a suitable sequence of RF and magnetic field 
gradient pulses. The result is an initial condition of NMR-labeled (i.e., polarized) ^^^Xe in 
the gas phase with unpolarized ^^^Xe in the dissolved phase. We then measure the diffusion 
of ^^^Xe polarization into the dissolved phase by monitoring the recovery of the NMR signal 
at the chemically-shifted frequency of dissolved-phase ^^^Xe. The rate of increase of the 
dissolved-phase ^^^Xe NMR signal is proportional to S/V (see analysis below). 

Fast, high-resolution NMR measurements of the nuclear spin-| noble gases ^He and 
^^^Xe are possible with optical pumping techniques using high-power lasers which increase 
the noble gas spin polarization to ~ 1 — 10% [|TB|. Laser-polarized noble gas NMR is a 



powerful technique in the physical and biomedical sciences, e.g., in lung imaging ||T^. For 
interphase exchange measurements, ^^^Xe is used rather than ^He because of two factors. 
First, ^^^Xe dissolves about ten times more readily than ^He in lipophilic compounds: e.g., 
^^^Xe has a large dissolved/gas partition coefficient of ~ 0.25 for lung tissue [|l^]. Second, the 
gas and dissolved-phase ^^^Xe NMR frequencies are typically separated (chemically-shifted) 
by ~ 200 ppm, a relatively large value (the typical '^He chemical shift is less than 1 ppm) 
which makes it easy to manipulate or observe ^^^Xe spins selectively in either phase ||19|| . 

Xenon interphase exchange NMR may be particularly important for the human lung 
where the mean pore (alveolar) size is ~ 100 /im 0,^. A noninvasive method to measure 
and image S/V and gas exchange in the human lung could be an important tool for probing 



lung physiology and pathophysiology pO[j ^.g., in the diagnosis and treatment of asthma 



and chronic obstructive pulmonary disease (COPD), the fourth leading cause of death in the 
U.S |21|]. Traditional lung function tests of gas exchange parameters are seriously limited 
by the inability to provoke step changes in gas composition or properties deep in the lung. 
Deconvolving the effects of ventilatory dilution, tissue absorption, and blood perfusion is a 
major challenge. By contrast, xenon interphase exchange NMR provides an almost ideal step 
change in ^^^Xe polarization between tissue and gas, and thus offers a unique window into 
lung structure and function. For example, subsequent to our annoucement of preliminary 
results using xenon interphase exchange NMR [|^, Ruppert et. al. succesfully applied a 
similar technqiue to studies of lung function in dogs p3| . Xenon interphase exchange NMR 
may also be a powerful tool in the characterization of non-living, soft porous media. For 
example, knowledge and control of S/V in packed columns of polymeric beads - commonly 
used in molecular fractionation applications - is critical to achieving adequate molecular 
separation. 

In the study reported here, we monitored the exchange of ^^^Xe polarization from the 
gas phase into the dissolved phase in a xenon-soluble porous polymer granulate (POREX 
Corp., Atlanta, GA), chosen as a model (with respect to mean pore size) for lung tissue. We 
investigated four different samples, which the manufacturer stated as having median pore 
sizes of 20, 70, 120, and 250 /im (obtained with mercury intrusion porosimetry on random 
production samples) and similar porosities of ~ 30%. The interphase exchange experimental 
protocol was as follows. First, we filled a glass container completely with a porous polymer 
sample, excluding free gas outside the cylindrical polymer core. Next, we placed the sample 
in an NMR magnet, pumped out the air from the sample's pore space, and allowed laser- 
polarized xenon gas to fill the space. (NMR imaging showed that the xenon gas uniformly 
filled the pore space within a few seconds.) We then used a frequency-selective RF pulse to 
quench the dissolved-phase ^^^Xe polarization, and observed the recovery of the dissolved- 
phase '^^^Xe NMR signal due to interphase exchange for times from 10 to 1000 ms. The NMR 
pulse sequence for the interphase exchange measurement is shown schematically in Fig. 
Note that we made an initial measure of the gas-phase ^^^Xe NMR signal using an RF pulse 



with a small flip angle a ~ 3 — 5 °. Also, the selective RF pulse at the dissolved-phase ^^^Xe 
NMR frequency was not perfect and hence generated a small NMR signal from ^^^Xe gas 
in the pores. We exploited these gas-phase signals to correct the data for slowly decreasing 
^^^Xe polarization due to both RF depletion and Ti relaxation in the pores. An example of 
the observed interphase exchange spectra is shown in Fig. |^. We performed xenon interphase 
exchange experiments on two NMR spectrometers at different magnetic field strengths: a 
horizontal-bore GE Omega/CSI spectrometer operating at 4.7 T (55.3 MHz ^^^Xe NMR 
frequency) with field gradient strength up to 7 G/cm; and a modified SMIS spectrometer 
interfaced to a prototype IBM horizontal-bore mag net operating at 1.5 T (17.6 MHz ^^^Xe 
NMR frequency) with field gradient strength up to 2 G/cm. For each experimental run, 
we used the spin-exchange optical pumping technique [|16| in the fringe field of the NMR 



magnet and achieved ~ 0.5% ^^^Xe spin polarization for 2 bar-liters of xenon gas (natural 
isotopic abundance). 

Calculation of absolute values for S/V required prior characterization of the samples. 
We used thermally-polarized xenon NMR on both spectrometers to determine the porosity 
and ^^^Xe partition coefficient b for the four porous polymer samples. We placed all 4 
sample cells on a manifold along with an empty cell of equal volume, admitted ~ 3 bar of 
xenon gas (natural isotopic abundance) and ~ 2 bar of oxygen, and allowed the system to 
equilibrate over 2-3 days to achieve full penetration of xenon into the polymer samples. We 
then measured the ^^^Xe NMR gas signal in the empty cell (s^o) as well as the equilibrium 
gas and dissolved-phase ^^^Xe NMR signals {sgas,Sdiss) in the four polymer samples. For 
each sample, = Sgas/sgo and b = {sdiss/ Sgas)(p/{^ — <P)- We obtained similar results for 
and 6 at 4.7 T and 1.5 T. Measurement of the dissolved-phase ^^^Xe diffusion coefficient 
Ddiss required large magnetic field gradients, which was only possible across small sample 
dimensions. Hence we cut additional pieces of each polymer sample to fit in standard 5 mm 
diameter NMR spectroscopy tubes and filled the tubes with ~ 10 bar xenon gas (isotopically- 
enriched to 90%^^^Xe). We then measured Ddiss using a 9.4 T Bruker AMX2 spectrometer 
with 50 G/cm field gradients, employing a modified pulsed-gradient technique previously 



applied to gas-phase diffusion measurements on rocks and glass beads [0 . Table | lists the 
measured values of 0, b, and Dctissi along with the ^^^Xe longitudinal spin relaxation times 
(Ti) for each porous polymer sample. 

To determine S/V from the xenon interphase exchange NMR measurements, we em- 
ployed a one- dimensional diffusion model, asymptotically valid for short times. We assumed 
that the density of dissolved-phase ^^^Xe polarization, pdiss{x,t), satisfies the diffusion equa- 
tion on a half space, x > 0, where the non-trivial spatial coordinate x is normal to the 
gas/solid interface. The quench of the dissolved-phase ^^^Xe polarization created the initial 
condition pdiss{x,0) = 0, x > 0. Given the much greater gas-phase ^^^Xe density, we also 
assumed a fixed boundary condition Pdissi^, t) = bpgas- Here, Pgas is the density of gas-phase 
^^^Xe polarization, which we found to be constant during an individual interphase exchange 
measurement except for the small effects of RF depletion and slow Ti relaxation mentioned 
above. The solution of the diffusion equation for Pdiss{x,t) is given by the well known error 
function |Q. The number of dissolved-phase polarized ^^^Xe atoms Ndiss{t) can then be 
obtained by direct integration of the density or by solving the ordinary differential equation 
obtained by integrating the diffusion equation. The result is: 

Ndissit) ,S 



y Pgasb—\/4:Ddisst/TT (l) 

where 5* is the total pore-space surface area of the sample, and the Euclidean volume V of 
the sample has been introduced as a common divisor on both sides of the equation. Noting 
that 

N as 

pgas = (2) 

where Ngas is the number of gas-phase polarized ^^^Xe atoms (assumed constant except for 
RF depletion and Ti relaxation), Eq. (1) can be rewritten as 

Ndissit) bS 



N,,. ^V4^-.V- (3) 

Ndiss{t) and Ngas are proportional to the frequency integrals over the observed ^^^Xe NMR 



spectral peaks, Sdissit,^) and Sgas{uj), for the dissolved and gas phases respectively. Thus 
we have 

J Sdiss(t,^^)duj bS 



J Sgas{uJ)duJ (f)V ^ 

We determined S/V by fitting the ratio of dissolved and gas-phase ^^^Xe NMR signals to 
\/t and using measured values for b, and D^iss 

Concurrent with the NMR studies, we also characterized the four samples of polymer 
granulate using confocal microscopy in order to give independent measurements of S/V. We 
carefully sliced sections of each of the four polymer rods with different pore sizes, creating 
glass smooth surfaces of both the longitudinal and cross- sectional faces of the rods. We 
oriented these sectional samples on a Leica TCS NT confocal microscopeconfigured for 488 



nm excitation and refiected light imaging [|25|. The resulting microscope images showed 



the polymer granulates as solids, and the pore spaces as black voids. Computer-based 
image analysis (ImageSpace, Molecular Dynamics, Sunnyvale, CA) superimposed multiple, 
equally-spaced lines across the images, and interface boundaries on each line of the images 
were counted. Thus for each sample we determined the mean linear intercept, L^, which 
is defined as the average distance between phase boundaries in an image. While is but 
one of several possible measures of pore size, it has the unique and powerful advantage of 
being stereologically related to the surface to volume ratio. In particular, is an unbiased 
estimator of S/V in the image plane: Lm = 2{S/V)^^ p6[. See Table H for the optically- 
determined values of Lm and S/V. 

We acquired three xenon interphase exchange data sets for each of the four different pore 
size polymer samples: two data sets at 4.7 T and one at 1.5 T. As an example. Fig. ^ shows 
an interphase exchange data set taken on the 4.7 T instrument. Each of the points plotted in 
Fig. ^ consists of the normalized average of three or four separate series of exchange spectra 
of the type shown in Fig. |^. The normalized xenon integral displayed on the ordinate of 
Fig. ^ is given by {I{t))(j)sina/b (see Eq. (4)), where a is the small flip angle of the RF 
pulse used in the initial gas-phase interrogation. All xenon interphase exchange NMR data 



followed the expected -y/t trend over the observed range of exchange times (out to ~ 1 s). 

Fig. ^compares S/V values derived from the optical and NMR techniques for each of the 
four polymer samples: one optical and three NMR data sets. (See also Table |l|.) The solid 
line in Fig. ^ is the line of identity for the optical and NMR values of S/V] the dashed line 
is a line of regression constrained to go through the origin. The slope of the regression line 
is 1.27 ± 0.18, indicating approximate agreement between the optical and NMR measures 
of S/V . Sources of systematic error in the interphase exchange NMR measurement include: 
miscalibration of the flip angle a of the gas-phase interrogation RF pulse; inaccuracy in the 
dissolved-phase xenon diffusion coefficient measurements; and deviations from the simple, 
one-dimensional diffusion model used to derive S/V . For example, the very low dissolved- 
phase xenon NMR signal from the nominally 120 /im POREX sample prevented us from 
measuring Ddiss for this sample. Instead, we used the average of the measured Ddiss values 
for the nominally 70 /im and 250 /im samples to calculate S/V from the interphase exchange 
measurements for the 120 /xm sample. 

We note that the manufacturer's quoted pore sizes for the four porous polymer samples 
do not agree with our pooled optical and NMR measurements of S/V (and hence L^); even 
the rank order does not agree. This discrepancy may not be surprising, since the pore sizes 
given by the manufacturer were determined by mercury intrusion porosimetry [P,^], which 
is primarily a probe of minimum pore-throat size and thus not necessarily a good measure 
of Lm or S/V p6[. We emphasize that there are different choices of the definition of pore 
size; which choice is most appropriate will depend on the particular application. 

In conclusion, we have demonstrated a new technique, xenon interphase exchange NMR, 
that provides a minimally invasive and nondestructive probe of the surface-area-to-volume 
ratio {S/V) of xenon-soluble porous materials. We have applied this technique to porous 
polymer samples and found approximate agreement with destructive stereological measure- 
ments of S/V made with optical confocal microscopy. 
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TABLES 

TABLE I. For each porous polymer sample, measured values are listed for the porosity 0, 
the ^■^''Xc gas/polymcr partition coefficient b, and the dissolved-phase ^^^Xe diffusion coefficient 
Ddiss- Also listed are the ^^^Xe longitudinal spin relaxation time (polarization lifetime) in the 
dissolved phase (Tiaiss) and the pore gas phase (Tigas)- The nominal pore size is that provided by 
the polymer manufacturer from mercury intrusion perimetry measurements on random production 
samples. Note: Ti^igg was the same at 4.7 T and 1.5 T; also, D^iss was not measured in the 
nominally 120 sample because the dissolved-phase ^^^Xe NMR signal from this sample was too 
low to allow the experiment to be completed in a reasonable time. Listed measurement uncertainties 
do not include estimates of systematic errors. 



Nominal pore size (/xm) 


20 


70 


120 


250 


Porosity, (p 


0.45 ± 0.02 


0.45 ± 0.03 


0.46 ± 0.01 


0.43 ± 0.02 


Xe partition coefficient, b 


0.64 ± 0.02 


0.52 ± 0.04 


0.56 ± 0.03 


0.63 ± 0.03 


XeDd^ss (x 10-12 m^s-i) 


5.6 ± 0.5 


7.0 ± 0.7 


n/a 


7.2 ± 0.7 


Xe Tidiss (s) at 4.7 T 


8.1 ± 0.4 


7.7 ± 0.3 


7.6 ± 0.3 


8.0 ± 0.3 


Xe Tigas (s) at 1.5 T 


21.8 ± 0.5 


34.1 ± 0.5 


35.4 ± 0.4 


19.9 ± 0.2 


Xe Tigas (s) at 4.7 T 


41.6 ± 0.6 


68.1 ± 1.4 


70.0 ± 1.6 


42.8 ± 0.6 



TABLE II. The surface-area-to-volume ratios (S/V) derived from the optical and NMR mea- 
surements are listed for each of the four polymer samples. Also listed is the mean linear intercept, 
Lfn = 2(5/^)-^, determined from the optical microscopy measurements. Listed measurement 

unccrtiiiulic^s do uol include; csiimalos of syslcmatic errors. 

Nominal pore size (//m) 20 70 120 250 

S/V from xenon NMR {nm'^) 0.030 ± 0.006 0.016 ± 0.005 0.016 ± 0.003 0.018 ± 0.003 
S/V from microscopy {nm'^) 0.0229 ± 0.0004 0.0135 ± 0.0003 0.0122 ± 0.0003 0.0127 ± 0.0003 
Lm from microscopy (/Ltm) 87.3 ± 1.7 147.8 ± 3.9 164.3 ± 4.5 158.0 ± 4.3 
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FIGURES 

FIG. 1. Xenon interphase exchange NMR pulse sequence. The sequence can be divided into 
the fohowing segments. A: interrogation of the gas-phase NMR signal {sgas) with a small flip 
angle RF pulse; Sgas is proportional to the ^^^Xe polarization. B: semi-selective saturation train 
of RF pulses to destroy the dissolved-phase ^^^Xe polarization. C: interphase exchange time, t. D: 
semi-selective 90° RF pulse to observe the recovery of the dissolved-phase ^^^Xe signal (sdissit))', 
this RF pulse also excites a small gas-phase signal, allowing correction for incremental loss of 
gas-phase magnetization. E: interphase exchange measurement, typically repeated ten times for 
each exchange time, t. 

FIG. 2. Dissolved-phase xenon NMR spectra as a function of interphase exchange time for the 
nominally 20 ^m pore-size sample. The peaks at ppm arise from the weak, far-off-resonance 
excitation of gas-phase polarized ^^^Xe when dissolved-phase ^^^Xe is pulsed semi-selectively. Fol- 
lowing saturation, the dissolved-phase ^^^Xe peaks at 200 ppm show recovery as a function of 
exchange time. The observed spectra are scaled to normalize the integrals of each gas peak to the 
gas peak in the first exchange spectrum; this normalization accounts for the slow loss of gas-phase 
magnetization due to exchange, RF depletion, and Ti relaxation. 

FIG. 3. One set of xenon interphase exchange NMR data, acquired at 4.7 T for each of the 
four porous polymer samples. (The samples are indicated in the figure by the manufacturer's 
claimed pore size.) Each data point represents the weighted mean of 3 to 4 normalized spectra 
of the type illustrated in Fig. ||, acquired in repeated runs of the experiment. The normalized 
xenon integral plotted on the ordinate is thus given by {I{t))(f)sina/b (see Eq. (4)). Uncertainties 
for each interphase exchange spectrum were determined from the measured rms noise within the 
spectral bandwidth, normalized for variations in gas-phase magnetization; one-standard-deviation 
error bars were then calculated for the weighted means of the normalized xenon integrals. The 
error bars are larger for longer exchange times because these data were acquired later in each 
experimental run and hence had lower gas-phase magnetization, which correspondingly increased 
the normalized noise of the dissolved-phase xenon spectra. 
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FIG. 4. Comparison of polymer pore-space S/V values obtained with optical microscopy and 
NMR. For each of the four polymer samples, three NMR experimental runs are plotted versus one 
optical measurement. The solid line is a line of unity for the S/V values provided by the two 
measurement techniques. The dashed line indicates a regression between the NMR and optical 
S/V values, constrained to pass through the origin; this regression line has a slope of 1.27 it 0.18. 
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